Fermion condensation around a Coulomb impurity in a Weyl semimetal and in a 
narrow band gap semiconductor as manifestations of the Landau zero-charge problem 
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A Coulomb impurity placed in an undoped Weyl semimetal spontaneously surrounds itself with 
a cloud of condensed Weyl fermions. We find that the ground-state of this system exhibits an 
experimentally accessible Landau zero-charge effect: the fermion condensate completely screens out 
the impurity charge. In a narrow band gap semiconducor this effect manifests itself in the near 
universality of observable charge of a highly-charged recombination center. 
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The Dirac equation for an electron in the field of 
a point charge Ze in vacuum becomes meaningless for 
Z > l/a « 137 (a = e^/hc is the fine structure con- 
stant) because it predicts an imaginary ground-state en- 
ergy [H- For a nucleus with Zc ~ 170 (upshifted from 
Z = 137 as it is not a point particle) the total energy 
of the production of an electron-positron pair vanishes 
and the vacuum becomes unstable with respect to pair 
creation; the positron repelled by the nucleus escapes to 
infinity while the electron remains near the nucleus 
For Z > Zf. vacuum electrons condense around the nu- 
cleus within a shell partially screening the nuclear charge. 
As a result, the observable dimensionless nuclear charge 
as seen at large distances is less than its bare value of Z 

Q • The peculiarity of this system is that its ground 
state, the vacuum, is charged. 

Here we observe that the prediction of vacuum elec- 
tron condensation in quantum electrodynamics (QED) 
and of "charged vacuum" can be tested in performable 
experiments involving condensed matter systems, both 
presently available and those that will become available 
in the near future. Our primary example is that of semi- 
conductors whose physics is known to mimic, to some 
extent, QED [2,]. The excitation of an electron- hole pair 
parallels the creation of an electron-positron pair in QED, 
with the band gap imitating the combined rest energy of 
the particles. There is also a counterpart to a .Z > 137 
nucleus in condensed matter physics. In his study of the 
impurity states in semiconductors Keldysh Q noted that 
the effective mass approximation ,6], while successful in 
describing shallow impurity states, fails to explain deep 
states whose binding energy is comparable with the band 
gap. Such states are formed near highly-charged impuri- 
ties and (in contrast to their shallow counterparts) they 
cannot be associated with either conduction or valence 
bands. The experiment presented another puzzle: some 
highly charged impurities acted as very efficient recombi- 
nation centers but an explanation why that was the case 
was lacking. Keldysh showed that experimental findings 
can be understood in a two-band approximation (well 
obeyed in narrow band gap semiconductors (NBGS) of 



the InSb type) where the low-energy electron-hole dis- 
persion law is relativistic [1, 0| 



e(p)=±y(A72)^ 



y2p2 



(1) 



where A is the energy band gap and v is the velocity of a 
high-momentum particle analogous to the speed of light 
in vacuum c; in NBGS we have v w 4.3 x 10~^c [8|. Then 
determination of the impurity states reduces to solving 
the Dirac equation in the field of a charge Ze screened 
by the dielectric constant e of the semiconductor. In 
view of the Za > 1 peculiarity of the Dirac equation 
(now a — e'^/hve), Keldysh concluded that for Za < 1 
the impurity states are given by the known solution to 
the Dirac equation jH while the anomalous Za > 1 case 
describes a recombination center. 

The recombination center is clearly a semiconductor 
counterpart of the Z > 137 nucleus. However, the semi- 
conductor equivalent of the fine structure constant is 
a = e^/hve ~ 1.1 /e. Since we now have e ~ 10 Q, 
the corresponding a is an order of magnitude larger than 
its QED counterpart. Thus the characteristic Z = l/a of 
the onset of the electron condensation is about 10. Surely 
a Z > 10 recombination center is a more common object 
than a > 137 nucleus. 

In addition to making it possible to study the regime 
of large effective fine structure constant, condensed mat- 
ter systems also offer possibilities that cannot be real- 
ized in QED. Indeed, over forty years ago Abrikosov and 
Beneslavskii [l3| predicted the existence of semimetals 
having points in the Brillouin zone where the valence 
and conduction bands meet with a dispersion law that 
is linear in the momentum. This is the A = case of 
Eq.(IT]). Such systems, nowadays called Weyl semimetals 
(WS), are likely to be realized in doped silver chalco- 
genides Ag^+gSe and Ag2+sTe [lH, pyrochlore iridates 
^2 ^''207 (where A is Yttrium or a lanthanide) Il2| , and 
in topological insulator multilayer structures [l3|. The 
zero energy gap of a WS implies that there is no thresh- 
old to creation of the electron-hole pairs. Thus conden- 
sation of Weyl fermions around a Coulomb impurity is 
spontaneous. The physically related problem of impu- 
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rity screening in graphene has been considered elsewhere 
|14| (that problem is mathematically different from what 
we discuss, because it is a two-dimensional semimetal 
embedded in a three-dimensional space). 

Below we will determine the ground-state properties 
of NBGS (including its WS limit) in the presence of a 
Coulomb impurity, as a function of Z and a. At modest 
Z the electrons of the "vacuum" (valence band) condense 
around the impurity while the holes leave the physical 
picture; the properties of the electron cloud vary with Z 
and a and are determined by the interplay of attraction 
to the impurity (promoting electron condensation), and 
electron-electron repulsion combined with the Pauli prin- 
ciple (limiting the condensation). The QED analysis of 
the physical properties of the vacuum electron cloud was 
carried out in two limits: 

(i) Z close to Zf.^ where there are very few condensed 
electrons and the single-particle picture is a good starting 
point 0, and 

(ii) Z 3> ^c, where the number of condensed electrons 
is large and the electron-electron interactions cannot be 
ignored 

The goal of this paper is to demonstrate that the 
physics of the charged vacuum in the Z Z^ limit ex- 
hibits large degree of universality manifesting itself in a 
nearly-universal observable impurity charge; its detection 
is within experimental capabilities. Although below we 
adopt the semiconductor language, our findings are also 
applicable in the QED implementation of the problem. 
Specifically, our central conclusion that in the gapless 
A = case the large distance character of screening is 
formally identical to that occurring in the Landau "zero- 
charge" problem [l[ was missed in previous studies Q. 

Since for Z Z^ & large number of condensed elec- 
trons is present, the properties of the system consisting of 
the impurity and its interacting cloud of electrons can be 
understood semiclassically with the help of the Thomas- 
Fermi (TF) theory d Q which becomes asymptotically 
exact in the Z oo limit 15|. The main object of the 
TF theory is a physical electrostatic potential (^(r) felt 
by an electron that is due to both the electrostatic po- 
tential of the impurity (^e^jt (r) and that of the condensed 
electrons characterized by the number density n(r): 



(^(r) = ^ext(r) - - / 



(2) 



The external potential ^Pextij) is a pseudopotential that 
represents the perturbation of the system caused by the 
impurity; even though kp^xt is not entirely of electrostatic 
origin, we will define l^i^ext = — 47renea;t/e. We assume 
that the impurity charge density euexti,'^) is spherically- 
symmetric and localized within a region of size a ~ Inm 
so that for r > a the potential ipext (r) reduces to a purely 
Coulomb form ipext [r) — Ze/er oi a net charge Ze within 
the impurity region. Given (p{r), one can deduce that 
the electron number density n{r) is different from zero 



only in the region of space where the electron potential 
energy —e(p{r) + A/2 drops below —A/2, thus defining 
the "vacuum" shell where condensed electrons reside as 



eip{r) > A, n(r) > 



(3) 



The radius of the electron shell Ryac is given by the equal- 
ities eip{Rvac) — A, n{Ryac) — 0; outside the shell we 
have n = and 

QooS n - "*5oo ^ eAa 

•■P = , r > R^ac = Qoo— r = Zq = (4) 

er eA Zq 

where Qoo < is the dimensionless observable impurity 
charge as seen at large distances from its center. Conti- 
nuity of the potential across the shell boundary relates 
Rvac and Qoo as indicated in the last two steps in (|4]) 
meaning that we can speak of the shell size or the ob- 
servable charge interchangeably. The parameter Zq < Zc 
gives the classical, a = oo, value of the critical charge of 
the onset of the electron condensation at the edge of the 
impurity region; in NBGS with A ~ O.leV and a ~ Inm 
we have ~ 1. In WS (Zq = 0) we find Rvac = oo, i.e. 
the screening cloud has infinite extent. 

From the thermodynamical standpoint, electron (e) 
condensation in the field of a Coulomb impurity accom- 
panied by escape of a hole (h) to infinity may be viewed 
as a " chemical reaction" e + h ^ (the ground state of 
the semiconductor is the "vacuum") [^; the condition of 
equilibrium for this reaction has the form 



0,/ie- v^(A/2)2 



eip,Hh = A/2 (5) 



where fie and fih are 
the electrons and holes. 



the chemical potentials of 
respectively, and PF(r) = 



h[6TT-^n{r) / g]^^^ is the Fermi momentum which we as- 
sume is a slowly varying function of position. The param- 
eter g is the fermion degeneracy factor: g = 2 in QED 
while in NBGS it is twice the number of Dirac valleys 
([T]) within the first Brillouin zone; an isotropic valley- 
independent limiting velocity v is assumed for simplicity. 
In the WS case g counts the number of Weyl points within 



the first Brillouin zone: (7 = 24 in pyrochlore iridates 12 1 
and g = 2 in a topological insulator multilayer [3] . Then 
the condition ([5]) implies a relationship between thephys- 
ical potential and the electron number density 0, |j] : 



n(r) 



A f eLp{r) e 



47r 



[eip{r) - A] 



3/2 



A = 



(6) 



where A is a dimensionless parameter that characterizes 
the relative strength of electron-electron interactions and 
zero-point motion. It can be as small as 1.65 x 10~^ 
(QED) or as "large" as lO"^ (NBGS). 

Applying the Laplacian operator to both sides of 
Eq.(l2]) and using ([6]) we find the relativistic TF equation 
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Figure 1: (Color online) Potential ([2| and charge Q[r) 
within a sphere of radius r ([13} as functions of distance (dou- 
ble logarithmic representation). The source region is r/a < 1. 
The gap value A is indicated by the horizontal line; the elec- 
tron cloud is limited to the region where > A, which de- 
fines Rvac- The curves are drawn for Z — 2, 20, 200, 2000, 
A = 0.001, and Zq ~ 1 (other values of Zo are equivalent to 
a rescaling Z — > Z/Zq and A — XZq). For large Z, Rvac 
approaches a Z-independent limit, indicating that Q{r) tends 
to an upper bound Qoo- 

that was investigated in QED 01 . The source term is 
locahzed and its specific form is not very important; we 
can take n^xt to have the constant value ^Z/Ana^. A nu- 
merical solution to ([7]) is shown in the Figure, where we 
additionally displayed the charge Q{r) within a sphere 
of radius r as an indicator of the strength of screening. 
In agreement with previous analysis [3, 01 the screen- 
ing effect of condensed electrons becomes noticeable for 
XZ^ > 1. Remarkably, for Z ^ 1 there exists a Z- 
independent limit on Ryac and Qoo- To understand the 
form of numerical solution in the strong screening regime 
AZ^ > 1, it is useful to start with the WS case, A = 0, 
when Eq.([7]) simplifies to 



(?) 



= -iTTUext + A 

V e 



(8) 



In the strong screening XZ^ — >■ 00 case the solution to 
Eq. dS]) is given by the zero of its right-hand side: 



n = no = next, 



tp = ipo 



1/3 



(9) 



We see that to leading XZ^ ~ 00 order the screening is 
complete with zero electric field everywhere and constant 
potential ipo inside the impurity region 0, 0] • In order 
to further improve on (jH]) , within the impurity region we 
substitute (p ~ (^o(l — /), < / ^ 1, into Eq.([9]) and 
linearize about ip ~ ipQ: 

K^f = 0, K^a^ = 3^/3(AZ2)i/3 > 1 (10) 



This approximation parallels the TF theory of screening 
in a Fermi gas [91]. We observe that inside the impurity 
region the screening response is characterized by the TF 
screening length ~ a{XZ^)~^^^ a: it is a length 
scale over which / drops to practically zero from a value it 
assumes on the impurity boundary. This is also the width 
of the region adjacent to the impurity boundary where 
uncompensated charge is localized [J|. From here the 
net charge inside the impurity region can be estimated 
as (Z/a3)a(AZ2)-i/6a2 ^ Z(AZ2)-i/6 < Z 0]. The 

crossing of the charge curves for Z = 200 and Z = 2000 
at small r shown in Figure is a direct consequence of 
screening: the TF screening length is smaller for 
Z — 2000 than Z — 200, so that the screening at the 
central region is more complete in the former case. 

The solution to Eq.([TU|) that is finite at the origin has 
the form / oc sinh(Kr)/r, so that for r ^ a. 



e e 



/) 



9Z 



I- A 



a sinh(Kr) 
r sinh(Ka) 



(11) 



where A is a constant assumed to be much smaller than 
unity to justify linearization approximation (jlOp . 

Outside of the source one has to look at the full non- 
linear equation whose solution is sought in the form 



r Va/ 



(12) 



where, via Gauss's theorem, the function C, is related to 
the charge Q{r) within a sphere of radius r as: 



Qir) = -r^^ 



C(0-C'(0, ^=ln- (13) 



Substituting into dS]) for r > a we obtain the equa- 
tion 



C"(0-C'(0-AC^ 



(14) 



For I — ln(r/a) ^ 1 we can neglect in Eq. p^ the 
first-order derivative term C'(Z) compared to C"(0; then 
Q{r) « -C'(0- The solution to dUl) in this limit is 



Ci(0 



1 



Xl + B' 



^ / < 1 



(15) 



Continuity of the potential and of the electric field at 
the impurity boundary r = a determines the integration 
constants A and B in Egs. pTj) and (|15l) to be 



0.2374, B 



1 



0.3113Ka 



(AZ2)-l/6 ^ ^ (;^g) 



The solution dlT]), ([T2|) and ^ also describes the NBGS 
case since the condition e(y5 3> A necessary for transition 
from ([7]) to ([8]) holds. Specifically, in the Z ^ 00 limit 
the parameter B vanishes and the Z-dependence drops 
out of din]). The solution to the full TF equation dl]) 
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for r > a then satisfies the singular boundary condition 
eip(r)/e — \/2/ X{r—a)~^ as r — )■ a leading to the numer- 
ically observed independent limit on i?„ac and Qoc- 

In the strong-screening limit the solution to the full 
Eq.dlll) has the form C(A, I) = {2/\f/'^y{l) where y{l) is a 
parameter free universal function such asy{l — > 0) — > l^^ . 
The latter behavior is no longer an accurate representa- 
tion of the true dependence y{l) past I ~ 1. Therefore 
the solution (fT5|) is only applicable up to a crossover scale 
Z = P ~ 1, i.e. within several impurity radii. Within this 
range the rescaled potential e(^/e drops from a value of 
the order A~^/'^(AZ^)^/^a~^ at the impurity boundary 
to X~^/''^a~^ at the crossover scale I* . This explains the 
large slope of the potential and charge curves near the 
boundary that can be observed in the Figure. 

For I — ln(r/a) 1 we can neglect in Eq. p^ the 
second-order derivative term C"(0 compared to C'(0; 
then Q{r) « and for arbitrary screening strength 
Eq. (fT4|) acquires a form 



(17) 



that is mathematically identical to the Gell-Mann-Low 
equation [l[ for the physical charge in QED reflecting 
the effects of vacuum polarization. Eq. ([T71) exhibits the 
effect of "zero charge" : no matter what the "initial" value 
of Q is, the system "flows" to the zero charge fixed point 
(5=:0as/— >oo,i.e. the impurity charge has been com- 
pletely screened. In the strong-screening regime Eq.(fT7|) 
is applicable at Z > Z* ~ 1. As a result the charge inside 
a sphere of radius r > a* — ae' > a will be given by 



1*2 



l + 2Ag*2in(r/a*) 2Aln(r/a*) 



(18) 



where the integration constant Q* is the charge within a 
sphere of radius a*. Since Q* ~ Ci(^*) — A"^/^ > 1, the 
charge Q{r) is given by the last representation in ([18]) 
whose hallmark is near universality: a weak logarithmic 
dependence on the source size a ~ a* with universal am- 
plitude. We conclude that in the WS case, except for 
the immediate vicinity of the impurity boundary where 
Eg. (lisp applies, the solution to the screening problem is 
nearly-universal. 

In the NBGS case (A ^ 0) the solution ([TH]) remains 
relevant at distances r <^ Rvac] however the true solution 
Q{r) corresponding to (fT5|) merges smoothly with Qoo 
at R vac- The latter can be estimated with logarithmic 
accuracy by equating eip to A, with the results 



1 



V2Alnr 



-^vac 



T 



r = 



1 



ZoV2A 



(19) 



The condensed electron cloud is of finite spatial extent 
and the screening is incomplete; the results (fTOl) are ac- 
curate provided InF ^ 1. For NBGS with a ~ Inm 
impurity region (Zg ~ 1) and A = 10^^ we find InF sa 3. 



This is not very large but sufficient to estimate the lim- 
iting charge as Qoo — (2A)~^/2 ss 22, and the size of the 
vacuum shell Rvac — 22nm, in agreement with numeri- 
cal solution of the problem. In the QED context we find 
Qoo ^ 137^/2 « 1600. 

We note that making the substitution Q* — > Z, a* ^ a 
in Eq. (fT8|) gives an accurate interpolation formula that 
describes the regime of weak screening {XZ'^ <C 1) for all 
r. The nearly universal limit of Eg. dTSl) will be reached 
at distances r > aexp{l/2XZ^). Applying to this so- 
lution the condition eip — A defining the edge of the 
electron vacuum shell, we recover, with logarithmic ac- 
curacy, earlier results Q corresponding to the regime of 
weak screening in NBGS. 

This work was supported by US AFOSR Grant No. 
FA9550-11-1-0297. 



[1] 

[2] 
[3] 
[4] 

[5] 



[6] 

[7] 
[8] 
[9] 



[10] 

[11] 
[12] 

[13] 
[14] 



V B Berestetskii, L. P. Pitaevskii, and E.M. Lif- 
shitz. Quantum Electrodynamics, Second Edition: Vol- 
ume 4 (Course of Theoretical Physics) (Butterworth- 
Heinemann, 1982), Sections 36 (QED Coulomb problem) 
and 133 (the Landau "zero charge" problem). 
Ya. B. Zel'dovich and V. S. Popov, Usp. Fiz. Nauk 105, 
403 (1971) [Sov. Phys. Uspekhi 14, 673 (1972)]. 

B. Miiller and J. Rafelski, Phys. Rev. Lett. 34, 349 
(1975). 

A. B. Migdal, D. N. Voskresenskii, and V. S. Popov, 

Pis'ma Zh. Eksp. Teor. Fiz. 24, 186 (1976) [Sov. Phys. 

JETP Lett. 24, 163 (1976)]; A. B. Migdal, V. S. Popov, 

and D. N. Voskresenskii, Zh. Eksp. Teor. Fiz. 72, 834 

(1977) [Sov. Phys. JETP 45, 436 (1977)]. 

L. V. Keldysh, Zh. Eksp. Teor. Fiz. 45, 364 (1963) [Sov. 

Phys. JETP 18, 253 (1964)]; see also P. A. Wolff, J. 

Phys. Chem. Sol. 25, 1057 (1964); A. G. Aronov and G. 

E. Pikus, Zh. Eksp. Teor. Fiz. 51, 281 (1966) [Sov. Phys. 

JETP 24, 188 (1967)], and M. H. Weiler, W. Zawadzki, 

and B. Lax, Phys. Rev. 163, 733 (1967). 

E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, 

vol.IX, Part 2, (Pergamon, 1980), Chapter VI. 

E. O. Kane, J. Phys. Chem. Sol. 1, 249 (1957). 

W. Zawadzki, Phys. Rev. B 72, 085217 (2005). 

C. Kittel, Introduction to Solid State Physics, 7th Edition 
(John Wiley & Sons, Inc., New York, 1996), Chapters 8 
and 10. 

A. A. Abrikosov and S. D. Beneslavskii, Zh. Eksp. Teor. 

Fiz. 59, 1280 (1970) [Sov. Phys. JETP 32, 699 (1971)]; 

J. Low Temp. Phys. 5, 141 (1971). 

A. A. Abrikosov, Phys. Rev. B 58, 2788 (1998). 

X. Wan, A.M. Turner, A. Vishwanath and S.Y. Savrasov, 

Phys. Rev. B 83, 205101 (2011); V. Aji, arXiv:1108.4426j 

A. A. Burkov and L. Balents, Phys. Rev. Lett. 107, 

127205 (2011); A. A. Burkov, M.D. Hook, L. Balents, 

Phys. Rev. B 84, 235126 (2011). 

D. P. DiVincenzo and E. J. Mele, Phys. Rev. B 29, 1685 
(1984); M. I. Katsnelson, Phys. Rev. B 74, 201401 (R) 
(2006); A. V. Shytov, M. I. Katsnelson, and L. S. Levitov, 
Phys. Rev. Lett. 99, 236801 (2007); M. M. Fogler, D. S. 
Novikov, and B. I. Shklovskii, Phys. Rev. B 76, 233402 



(2007). 

[15] L. Spruch, Rev. Mod. Phys. 63, 151 (1991). 



